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ABSTRACT 



in 



Unidentified extreme UV sources, detected in the EUVE and ROSAT WFC all-sky 
surveys, could be isolated old neutron stars, accreting material from the ISM. The 
closest neutron stars, that are located in the local ISM bubble of unusually low density, 



faint and cool (L 



10 



27 



T < 6eV). The extreme UV spectrum of these 
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O About 10% of the sources detected in the all sky surveys of 
£jthe EUVE (Malina et al. 1994, Bowyer et al. 1994) and the 
ROSAT WFC (Pound et al. 1993) have no previous ident- 
ification at optical or other energies. The nature of these 
sources, generically classified as 'NOIDs', is unclear and 
probably inhomogeneous. There are no strong constraints 
on their distance, some of the NOIDs may be extragalactic 
and some may contain variable and flaring sources. However, 
there is the exciting possibility that some of them might be 
nearby isolated old neutron stars (i.e. neutron stars not in 
binaries or in supernova remnants, and too old to emit ra- 
diation as normal pulsars) accreting material from the ISM 
(Ostriker, Rees & Silk 1970, Treves & Colpi 1991, Blaes & 
Madau 1993). Our Galaxy is estimated to have ~ 10 9 iso- 
lated old neutron stars, a population 3-4 orders of magni- 
tude larger than the pulsar population, but that has so far 
received little attention. The closest isolated neutron star 
is likely to be less than ~ 10 pc from Earth, and several 
hundreds are expected within 100 pc. 

The EUVE survey, lasted from 24 July 1992 to 21 Jan- 
uary 1993, covers 97% of the sky with exposures ranging 
from a few hundred seconds at the ecliptic to 20 kilosec- 
onds at the poles. The survey was carried out in paral- 
lel with three scanners, and a Deep Survey/Spectrometer 



sources is very sensitive to the H I column density, since large fraction of the emission 
is emitted just below the hydrogen Lyman edge. The EUVE sources with large count 
rate in the long-wave bandpass (600A) seem to be the most promising candidates. 
These sou 



should have low H I column density (Nhi ^ 10 cm ), constraining 
their distance to a few tens of parsec. Otherwise, their spectra would be significantly 
modified by ISM absorption, inevitably become stronger in the short wavelength (100A, 
200A) EUVE bandpasses. If these unidentified objects are familiar EUV sources rather 
than neutron stars, either white dwarfs, late type stars or cataclysmic variables, they 
are expected to be identifiable, generally brighter than V ~ 14. 
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(DS/S) telescope, pointed in perpendicular to the scanners. 
The peak wavelength and wavelength regions at 10% filter 
bandpasses are 

1) Scanner A & B, filter Lexan/B ("Lex") A = 891 
AA = 58 — 1741; filter Al/Ti/C ("A1C") A = 17ll AA = 
156 - 2341; 

2) Scanner C, filter Al/Ti/C/Sb ("Dagwood") A = 
4051 AA = 345 - 6051; filter Si/SiO ("Tin") A55ll, 
AA = 519 - 7431; 

3) DS/S, filter Lexan/B A = 9ll AA = 67 - 1781; filter 
Al/C A = 17ll, AA = 157 - 3641. 

The ROSAT survey, started on July 1991 and lasted 
about half a year, was performed with two filters, SI, AA = 
60 - 1401, and S2, AA = 110 - 200l. 

The unidentified EUV sources detected by EUVE and 
R OS A T WFC are generally faint, but rather inhomogeneous 
in terms of count rates in each wavelength range. The spa- 
tial distribution of NOID sources is only weakly concen- 
trated toward the Galactic disk. In this study we concentrate 
on a subclass of EUVE NOID sources, whose flux in the 
longest wavelength band (Tin) is remarkably strong. Spe- 
cific sources associated with this class are EUVE J0256+080, 
EUVE J1244-596, EUVE J1533+337, EUVE J1636-285 and 
EUVE J1748+484 (The 2nd EUVE catalog). The hydrogen 
column density estimated for the space positions of these 
sources suggest that none of them is likely to be more than 
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lOOpc away (B. Welsh, private communication). For neutral 
hydrogen column density Nhi = 10 17 cm -2 the ISM trans- 
mission in the Tin band is ~ 0.7, the same order of magni- 
tude as the transmission in the Lex and A1C bands (~ 1). 
However, when the column density is raised to 10 18 cm -2 , 
the Tin-to-Lex transmission ratio is lowered to ~ 0.1, and 
when Nhi = 10 19 cm -2 this transmission ratio reduces to 
~ 10 -10 (!). Therefore, independent of their spectral distri- 
bution, the column density towards sources of 'Tin excess' 
flux is strongly constrained, to Nhi ~ 10 18 cm -2 . 



2 THE EUV COUNT RATES OF AN 
ACCRETING NEUTRON STAR 

In hydrodynamical spherical accretion (Bondi & Hoyle 1944) 
the accretion rate is strongly dependent on the inverse of 
the star's velocity and the density of the ambient mat- 
ter (M cx nv~ 3 ). The average speed of Galactic neutron 
stars (Madau & Blaes 1994) is believed to be rather large, 
~ 100 km s -1 . Only matter enclosed within the accretion ra- 
dius face ~ 3.7 x 10 12 «^o cm (where v\oo = v/100 km s -1 ) 
could be accreted onto the neutron star. The neutral gas 
density in the local interstellar medium is unusually low, 
estimated to be < 0.01cm -3 for most directions out to 50 
pc (Cox & Reynolds 1987, Welsh et al. 1994). Accretion 
onto close neutron stars would therefore occur on low rate. 
Moreover, the accretion efficiency is even further lowered if 
the particles are not coupled, where the fluid approximation 
breaks down. In our case the gas is indeed collisionless , since 
face is orders of magnitude shorter than mean free path for 
Coulomb collisions (~ lO 18 ™,^ cm; rao.oi = ra/0.01 cm -3 ). 
However, the particles, that are photoionized by emission 
from the accreting neutron star, are coupled by ambient 
magnetic fields and plasma instabilities (for B = 10 -6 G 
the Larmor radius is 3-4 orders of magnitude smaller than 
the accretion radius). 
Low accretion rate 

M = 2.9 x 10 r ra .oi«roo g s_1 (!) 

could yield bolometric luminosity of only a small fraction of 

Lq, 

■ GM „ ■ , , 

L = f em M = f em 1.8 x 10 27 M 7 r 6 -1 erg s -1 . (2) 

r 

Here f em is the fraction of the accretion luminosity being 
converted to radiation, M is the neutron star mass, taken 
to be 1.4Mq, re is the radius in units of 10 6 cm and Mr = 
M/10 7 g s -1 . 

The accretion power can be converted to radiation 
through bremsstrahlung (mostly ion-electron, as T < 0.51 
MeV), Compton and, in the case of strong magnetic fields, 
synchrotron processes. Treves & Colpi (1991), Blaes & 
Madau (1993), and Madau & Blaes (1994) assume the ac- 
creted matter is channeled by the dipole magnetic field to- 
ward the magnetic poles and the emission is released from 
optically thick plasma in the star's atmosphere above the 
polar caps. Nelson et al. (1994) argue that a few percent of 
the luminosity in that case should be nonthermal, emitted in 
Comptonized cyclotron photons just below 11.6(B/10 12 G) 
KeV. In an alternate, optically thin model, particles are ac- 
celerated at the Alfven surface, far from the neutron star 



in a standing shock wave. The energetic particles upscatter 
the ambient soft photons and emit synchrotron radiation. 
The resulting emission spectrum is nonthermal, generically 
a power law (Shemi 1995). 

The intrinsic spectrum of an EUV source is modified 
by the attenuation of the ISM and by the filtering of the 
telescope. Denoting F the intrinsic photon spectral distri- 
bution, TR the ISM transmission, and A the telescope ef- 
fective area, the resulting spectrum of the source would be 
F(X)TR(X)A(X). For a source at distance d the overall count 
rate in a bandpass AA is the integral over this spectrum, 

count rate = — [ — / F(X)FR(X)A(X)dX . (3) 
Alrd Jax 

In this study we have calculated the count rates in the 
four EUVE bandpasses, assuming the radiation is released 
from optically thick regions above the neutron star mag- 
netic polar caps. The blackbody effective temperature in 
that model is given by 

T eff ~ 20/ e 1 / I 4 M 7 1/4 A -1 / j t eV (4) 

where A caPl i is the emitting area, in units of 1 km 2 . 

Figure 1 show spectra of optically thick sources with 
blackbody at T e // = 2, 4, and 6 eV, respectively. For each 
temperature we have calculated several cases, corresponding 
to variable column density from to 10 19 cm -2 . The three 
cases in the left side of the figure are related to small H I 
column density Nhi = 10 17 cm -2 . These spectra found to 
be similar to spectra obtained with a vanishing Nhi, which 
means that the blackbody curve is mostly modified by the 
telescope filtering. The cases in the right side correspond to 
strong ISM attenuation, Nhi = 10 18 cm -2 . Table 1 show 
the wavelength-integrated count fluxes (equation 3) for the 
four Nhi values we have used (numerical values in Fig 1 
and Table 1 are normalized with a Tin count rate of 400 
count/Ksec). The data for effective areas vs. wavelength, 
and models of ISM transmission vs. wavelength and hydro- 
gen column density, have been taken from the EUVE cycle 
III NRA. The data is available for the public in the EUVE 
WWW node. 



3 DISCUSSION 

The most prominent conclusion arises from Figure 1 and 
Table 1 is that the 'Tin excess' NOIDs are likely to be 
rather cool, T e // < 6eV and have low column density, N < 
10 18 cm -2 . The distance of these sources, unless are located 
in an extremely rarefied 'tunnel' in the ISM, is constrained 
to d < 30 ra^Jj pc. Moreover, as is demonstrated in the cases 
of T = 6 eV, even for negligible ISM attenuation our model 
excludes hot sources (T ^ 6 eV). 

The bolometric magnitude of a source with L = 2 x 
10 27 erg s -1 is 19.75. Using T eff = 25,000 K° (Ft — V = 
— 0.2) and bolometric correction rriboi — m v = 2.6, we obtain 
the V (B) magnitude 22.26 (22.06) + 5 log d 10pc (dio = 
d/10 pc). The Tin count rate of such a source is found to be 
~ 50 d^g count/Ksec. Taking these values, we can predict 
the distance of EUVE J0256+080 (Tin=256 ± 57), EUVE 
J1244-596 (199 ± 43), EUVE J1533+337 (131 ± 34), EUVE 
J1636-285 (410 ± 79) and EUVE J1748+484 (133 ± 32) to 
be ~ 35, 40, 50, 28 and 50 x M T pc, respectively. 
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Figure 1. Spectra of accreting neutron stars, for Njjj = 10 17 
(left) and 10 18 cm -2 (right). Upper: T = 2 eV, central: T = 4 
eV, lower: T = 6 eV. Spectra corresponding to the Tin filter are 
plotted with solid lines. Spectra of Lex, A1C and Dag are drawn 
with dotted, short-dashed and dotted-dashed lines, respectively. 
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Table 1. Count rates in EUVE bandpasses 
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So far we have taken the conservative approach and con- 
sidered sources whose spectra are blackbody. Photoioniza- 
tion of He I and He II in the source atmosphere could, how- 
ever, reduce the photon flux below the He I and He II Lyman 
edges (504A and 228A, respectively) and significantly reduce 
the EUV-to-UV flux ratio. Such continuum 'jumps', of about 
two orders of magnitude, are predicted in models of early 
type stars (Kurucz 1979). The spectrum of the strongest 
EUV source detected so far, the B2 II star e CMa, indeed 
confirms atmospheric He I and He II attenuation, although 
the observed jumps are shallower than predicted (Cassinelli 
et al. 1995). In order to take into account such an atmo- 
sphere, whether viable or not, we included synthetic jumps 
in spectra that are inputed in our calculations. We found 
only a moderate effect in the spectrum. The ISM attenua- 
tion, once it is not negligible, is always the most important 
factor in modifying the blackbody curve. The conclusion is 
that, unlike accreting neutron stars in binary systems, hot 
(T e ff ~ 100 eV) sources remain improbable. 

Using equation 4 one can see how the low temperature 
we obtained strongly constrains the physical parameters of 
the source, f em Mr A~^ p 1 < 10 -2 . With this constrain, as- 
suming large efficiency f em ~ 1 and plausible accretion rate 
Mr > 1, the emitting area required to be much larger than 
in X-ray binaries. In X-ray binaries the emission is believed 
to be released from an accretion column, that is supported 
by a strong dipole magnetic field. If the emitting area in 
our case is hundreds time larger, the magnetic field should 
be significantly weaker. For T e // < 3 eV Eq. 4 translates 
to an emitting area comparable with a neutron star sur- 
face, > 10 13 cm 2 . One can infer that the flow configura- 
tion remains spherical even near the neutron star. Mag- 
netic field pressure is insufficiently strong to balance ram 
pressure of the flow when approaching the star. When the 
flow is converged the ram pressure is enhanced, roughly by 
(racc/r 6 f /2 ~ 3xl0 16 . From the pressure imbalance we find 
an upper limit on the magnetic field B < 750raJ / ' o 2 1 v t gj 2 G. 
Such a field is ten orders of magnitude (!) lower than mag- 
netic fields in pulsars. Evidently, we consider this value with 
great caution and leave it to be studied elsewhere, as the 
detailed physics of the accretion stands beyond the scope of 
this Letter. Nevertheless, it is worth remarking that such a 



low magnetic field, if real, could be evidence for strong field 
decay. 

If they are not neutron stars, what could the EUV 'Tin 
excess' NOIDs be? Assume the source of the EUV radiation 
is an early type star like e CMa, with similar count rate ra- 
tios and optical luminosity (e CMa : Lex : A1C : Dag : Tin 
= : 80 : 4030 : 41840, m v = 1.5, T eff ~ 21, 000 K°). Simple 
scaling gives a bright m, ~ 6.5 star. As no optical counter- 
parts brighter than m, = 14 appear within the error circle 
of the NOID sources mentioned above, early type stars are 
unlikely. Other EUV sources like white dwarfs or late type 
stars, if having M v < 14, are also improbable, while exten- 
sive optical observations of the NOID fields are required to 
constrain fainter optical counterparts (Maoz et al. 1995). 

Selection and a comprehensive study of best neutron 
star candidates among the NOID sources need deep optical 
UV, EUV and X-ray observations. Simultaneous observa- 
tions, for instance with the future Russian - International 
Spectrum-X- Gamma observatory, will be very useful to con- 
strain eruptive sources, like flare stars or cataclysmic vari- 
ables. Moreover, repeat of observations with good angular 
resolution, on time scale of a few years, could confirm proper 
motions of 2 d^v\oosin0 yr -1 , where 9 is the angle be- 
tween the line of sight and the star trajectory. 

The detection of isolated old neutron stars will provide 
extremely valuable data on the neutron star physics and 
evolution. Since isolated neutron stars are the fossils of the 
high-mass end of the initial mass function, their detection 
could provide important clues for deciphering the stellar evo- 
lution history of the Galaxy. Isolated neutron stars may also 
be important as Galactic gamma-ray burst sources (Higdon 
& Lingenfelter 1990), cosmic ray accelerators (Shemi 1995), 
and the origin of the Galactic excess in the X-ray background 
(Maoz & Grindlay 1994). 



Acknowledgments 

Astronomy at the Wise Observatory is supported by grants 
from the Ministry of Science and the Arts and from the Israel 
Academy of Science. The author acknowledges discussions 
with Piero Madau, Dan Maoz and Barry Welsh. 



REFERENCES 

Blaes, O. & Madau, P. 1993, ApJ 403, 690 

Bondi, H. & Hoyle, F. 1944, MNRAS 104, 273 

Bowyer, S. et al. 1994, ApJ Suppl. 93, 569 

Cassinelli, J. P. et al. 1995, ApJ, in press 

Cox., D. P. & Reynolds, R. J. 1987, ARA&A 25, 303 

Higdon, J. C. & Lingenfelter, R. E. 1990, ARA&A 28, 401 

Kurucz, R.L. 1979, ApJS, 40, 1 

Madau, P. & Blaes, O. 1994, ApJ 423, 748 

Malina, R. F. et al. 1994 AJ 107, 751 

Maoz, D. et al. , in preparation 

Maoz, E. & Grindlay, J. 1994, CFA Preprint. 

Nelson, R.W., Salpeter, E.E. Wasserman I. 1993, ApJ 418, 874 

Ostriker, J.P., Rees, M.J. & Silk, J. 1970, Ast. Lett. 6, 179 

Pounds, K. A. et al. 1993, MNRAS, 260, 77 

Shemi, A., 1995, MNRAS, in press 

Treves, A. & Colpi, M., 1991, AA 241, 107 

Welsh B. Y., Craig N., Vedder, P. W. & Vallerga, V. V 1994, 
ApJ, 437, 638 



